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Abstract 

Magnetic impurities coupled antiferromagnetically to a one-dimensional Heisenberg 
model are studied by numerical diagonalization of chains of finite clusters. By calcu- 
lating the binding energy and the correlation function, it is shown that a local singlet 
develops around each impurity. This holds true for systems with a single impurity, with 
two impurities, and for impurities forming a lattice. The local character of the singlet is 
found to be little affected by the presence of other impurity spins. A small effective inter- 
action is found between a pair of impurity spins, which oscillates depending on impurity 
distances. For impurity lattices, the energy spectrum shows a gap which is found to be 
much smaller than the binding energy per impurity if the coupling constants are small. 
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For larger coupling constants, it increases to the same order of magnitude as the binding 
energy, indicating that a local singlet is broken to create excited states. Impurity lattices 
with ferromagnetic couplings are also studied and their connection to the Haldane problem 
is discussed. 

PACS numbers: 75.10Jm, 75.20.Hr, 75.30.Hx, 75.30.Mb 
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I. INTRODUCTION 

A magnetic impurity embedded in a system of conduction electrons forms at zero tem- 
perature a spin singlet through the antiferromagnetic Kondo coupling. This Kondo prob- 
lem has been extensively studied and even solved exactly. ^'^ In these studies, the interac- 
tion between conduction electrons is usually neglected. The general question arises what 
changes we should expect when the interaction between conduction electrons becomes im- 
portant so that they are strongly correlated. As a first step toward an answer of this 
question, we consider strongly interacting conduction electrons in a half-filling band in one 
dimension. Charge fiuctuations are suppressed by the strong interaction, and the physics 
of the low-lying excitations is well described by the spin-| antiferromagnetic Heisenberg 
model in this strong coupling limit. The Hamiltonian is given by^ 

H = H, + Hfa, (1.1) 

where 

Hc = JY1 i^icS^c + o^{SlS% + SlSl)] , a = 1, (1.2) 

<i,j> 

Hfc = J'Sqc ■ So/, (1.3) 

with J > 0, J' > 0. Here Sjc and Sq/ denote the spin-| operators for the electron at site 
i of the chain and for an electron in an impurity which is close to site 0. The notation 
< i,j > refers to nearest-neighbor sites in the chain. A periodic boundary condition is 
imposed on the chain system. The system is schematically shown in Fig.l. 

In the presence of an Ising-like anisotropy (ct < 1 in Eq.(1.2)), the ground state is 
Neel ordered, and the excitations described by He have an energy gap. In the classical 
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picture, the impurity spin is antiparallel to the neighboring spin in the chain. Quantum 
fluctuations reduce | < > | from its classical value, but do not completely destroy the 
classical picture. [Here < A> denotes the average of A with respect to the ground state.] 
As analyzed in the Appendix on the basis of a spin-wave expansion,^'^ | < Sqj^ > \ decreases 
rapidly when a approaches 1, while the transverse spin correlations | < Sqj^Sq^ > \ grow 
to nearly maximum value 0.25. Although the spin- wave expansion may not work well for 
a ^ 1, these results suggest that in this limit the quantum fluctuations are large enough 
for the impurity spin to form a singlet with the spins of the chain. The purpose of this 
paper is to study the nature of this singlet formation for a = 1 by exact diagonalization 
of flnite-size clusters, which has been successful to study quantum spin models in one 
dimension.^ In order to obtain a comprehensive understanding of the problem, we consider 
not only the single impurity case but also the cases when two impurities are present or 
when the impurities form a lattice (see Fig.l). 

If one of the impurity spins forms a singlet with the spin in the chain closest to the 
impurity site, the energy gain by Hfc is | J'. This local singlet, however, will lose some 
energy due to the presence of He- Calculating the binding energy, we flnd that the energy 
gain per impurity spin is considerably reduced from As shown later, the reduced 

values are, however, still much larger than those for the conventional Kondo problem,^ 
suggesting that the binding energy is enhanced by increasing interactions between con- 
duction electrons. It is also found that the binding energy per impurity spin is nearly 
independent of the number of impurities. 

We also calculate the correlation function between an impurity spin and the spin in the 
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chain closest to the impurity site. The absolute value is found to increase rapidly with 
increasing J' . Just like the binding energy, it is nearly the same for the three different 
systems under consideration. Even for small values of J', this coincidence seems to hold. 
The local disturbance generated by an impurity site seems to be little influenced by the 
presence of other impurities. However, this does not necessarily mean that correlations 
between impurities are small; as shown later, they may be substantial, particularly for 
small values of J' . The reason is that a local triplet is admixed to the singlet state due 
to interaction with the other spins in the chain and this results in an indirect interaction 
between impurities and an increase in the binding energy. In the conventional Kondo 
model, the two-impurity problem has been studied by several methods. ^"^"^ Thereby a 
main issue is the interplay between the single-site Kondo effect and the RKKY interaction. 
In particular, the half-filled one-dimensional Kondo lattice has been studied by several 
authors. ■'^■'^"■'^^ It is known that for small Kondo-coupling constants the RKKY interaction 
becomes more important than the Kondo-singlet formation. The present findings extend 
those observations to the case of strongly correlated electrons. 

In the impurity-lattice case, we also calculate the energy gap between the lowest excited 
states and the ground state. It is found that they are triplets with momentum k = tt when 
the lattice constant is set equal to one. The gap is very small for small values of J' . This 
is consistent with a non-analytic dependence oc e~°'/'^' where a is some constant which is 
found for the half-filled one-dimensional Kondo lattice. The gap increases with increasing 
values of J' , and eventually becomes of the order of the binding energy per impurity spin. 
This suggests that the excited states are local triplets which are created by breaking up 
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local singlets. This result is consistent with the one found for the double chain model, ^ 
which is slightly different from the present model. 

In addition to the studies involving an antiferromagnetic coupling, we also study systems 
of impurity lattices with a ferromagnetic coupling (J' < 0) to the spins in the chain. By 
calculating the binding energy and the correlation function, it is shown that local triplet 
states develop with increasing values of |J'|. As a whole, the ground state is a singlet, 
and the lowest excited states are triplets with momentum k = tt, just as before. The 
energy gap is very small, and increases gradually with increasing values of | J'|. It is much 
smaller than the binding energy per impurity spin, indicating that the excited states do 
not destroy the local triplets which are mainly responsible for the binding energy. In the 
limit J' — > — oo, the gap goes over into the Haldane gap.-^^ In this context, there have 
been several models studied which are slightly different from the present one, such as the 
ferromagnetic-antiferromagnetic alternating Heisenberg chain model^^ or the double chain 
model. The present results are consistent with the ones found previously. 

In Sec. II, numerical analysis is given of exact diagonalization of finite clusters. Section 
III contains the concluding remarks. In the Appendix, a system is analyzed within the 
spin-wave expansion, in which an impurity spin is coupled to an Ising-like anisotropic 
chain. 

II. NUMERICAL ANALYSIS 

In the following we present the results for binding energies as well as various pair 
correlation functions of finite Heisenberg chain systems with a single impurity, with two 
impurities, and with impurities forming a lattice.^® Thereby we measure energies in units 
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of J. 



A. Single impurity 



We consider the Hamiltonian given by Eq.(l.l) with a = 1. The ground state is a 
singlet. The binding energy A£'i(J', N) is defined by 



where Eg{J',N) denotes the ground-state energy for a spin chain consisting of sites 
with an impurity coupled to it according to Hfc- The number of sites N varies from 4 
to 23. Values for AEi{J', N) are naturally grouped into two series according to even or 
odd numbers of N. The values for odd numbers N are larger than those for even numbers. 
This is explained as follows. The ground state for chains with odd numbers N is a doublet. 
From it a localized doublet may be easily generated to form a singlet with the impurity 
spin. But the ground state for chains of even numbers N is a singlet, from which a localized 
doublet cannot be created. With increasing N, binding energies decrease for N being odd, 
while they increase for even numbers of N. Since effects caused by the impurity spin are 
expected to be confined to a finite range from the impurity site, one might think that the 
binding energy is independent of the system size if it exceeds the size of the spin correlation. 
This is not the case though, since the energy gain due to singlet formation depends on the 
energy spectrum of the chain, which is discrete and depends on the system size. When 
extrapolating each series to A/" — > oo, we assume a dependence of AEi in powers of 1/N, 



AEi(J', N) = EgiO, N) - Eg{J', N), 



(2.1) 




(2.2) 
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to which we fit the data by a least-square method. Since the two series approach each 
other, we can evaluate accurately their limiting value. If the impurity spin forms a per- 
fect singlet with the spin in the chain closest to it, the energy gain due to Hfc would be 
I J'. This singlet state suffers an energy loss through He though. Therefore some com- 
promise between the two energies must take place. We find that the net binding energy is 
considerably reduced from | J', with a factor which is larger for smaller values of J'. 

We also calculate the correlation function < Sq/ • Sqc > between the impurity spin 
and the spin in the chain closest to the impurity site. Figure 3 shows calculated values 
of — < So/ • Soc > as function of J'. Like for the binding energy, they fall into two 
series according to odd and even values of N. The correlation is antiferromagnetic. With 
increasing N, — < So/-Soc > decreases for odd values of N, while it increases for even 
values, just like in the case of the binding energy. We extrapolate each series to A?" — > oo in a 
similar way as done for the binding energy, although the assumption of a power expansion 
in 1/N is uncertain for this quantity. The values extrapolated by the two series are so 
close that we obtain reliable estimates of the correlation function for J' > 0.2. For small 
values of J' (< 0.2), however, the two series do not approach each other closely enough 
to obtain reliable estimates. This suggests that the assumption concerning the power 
expansion in 1/N is inadequate. The absolute values first increase rapidly and then move 
gradually toward 0.75 with increasing values of J', which is the value for a perfect singlet. 
Contributions from triplet states of the two spins So/ and Soc appreciably infiuence the 
value of < So/ ■ Soc > provided J' is not too large. Another way of stating this is by saying 
that the impurity spin forms a singlet with the spins of the chain over a region which 
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is rather extended. Spin correlations between an impurity spin and that of conduction 
electrons in a Hubbard chain have also been considered by Hallberg and Balseiro,^^ who 
found an oscillatory behavior with distance. 

B. Two impurities 

In this subsection, we study the case of two impurities placed on nearest-neighbor and 
on next- nearest-neighbor sites, respectively, as shown in Fig. 1(b). 

The binding energies, AE2'^{J', N), AE2^^{J',N) are defined by the same equation 
as Eq.(2.1), where superscripts nn and nnn stand for nearest-neighbor and next- nearest- 
neighbor positions, respectively. Figure 4 shows the binding energy per impurity spin 
AE2^{J', N)/2. The series of values for even and odd numbers of sites N approach each 
other more rapidly than for the single-impurity case. We have also calculated 
AE^"^(J',iV)/2, whose dependence on J' and N is similar to AE^"(J', iV)/2. Figure 
5 shows the extrapolated values to — > oo for AE2^{J')/2, A£'2""(J')/2 and their com- 
parison with AEi{J'). Both values are close to A£'i(J'), indicating that the presence of 
another impurity spin affects the first one so little that the impurities can be considered as 
being almost independent. The value AE2^{J') /2 is slightly larger, while A£'2""(J')/2 is 
slightly smaller than AEi{J'). The mutual influence of impurities, though it is small, seems 
to show a small energy gain or loss respectively, i.e., an oscillatory behavior, depending on 
the distance between the two impurities. 

Figure 6 shows the correlation function — < So/ • Sqc > between one of the impurity 
spins and the spin in the chain closest to it, when the two impurities are placed onto 
nearest- neighbor sites. A local singlet develops with increasing values of J' . By comparing 
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Fig. 6 with Fig.3, we find that the extrapolated values are very close to those for the 
single-impurity case. Therefore the presence of second impurity has little effect on the pair 
correlation, in agreement with the findings for the binding energy. 

The small mutual effect on both, the binding energy and the correlation function < 
So/ • Soc >, suggests a small correlation < Sq/ • Si/ > between two impurity spins. Figure 
7 shows the calculated values of — < So/ • Si/ > for impurities placed on nearest- neighbor 
sites. Against what we would expect, their absolute values are not small. An impurity 
spin does not form a perfect singlet with the spin in the chain closest to it, instead a 
triplet state is mixed in, as shown in Fig.6. We speculate that an indirect interaction, 
though small, is based on these degrees of freedom, thereby affecting considerably the spin 
correlation. The correlations decrease monotonically with increasing values of J' for even 
numbers N, while they first increase and then decrease for odd numbers N. Their decrease 
with increasing J' (J' > 0.25) may be taken as an indication of the development of a 
singlet state between each impurity spin and the spins in the chain, resulting in weaker 
correlations between impurity spins. The values extrapolated to A?" ^ oo for even and 
odd values of N are sufficiently close in order to allow for reliable estimates of the limiting 
behavior provided J' > 0.25. But when J' — > 0, the values for even and odd numbers of N 
differ considerably. This may be interpreted as follows. In the absence of Hfc, the singlet 
and triplet states of two impurity spins are degenerate. Now let us turn on Hfc- If J' is 
smaller than the finite gap in the excitation energy of the finite-size systems, the state of 
the spins in the chain is little modified by the impurity spins. For even values of N, the 
ground state of the spins in the chain is 'nearly' a singlet. Since the total system must be a 
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singlet, the singlet state between two impurity spins is favored, resulting in large absolute 
values of the pair correlation. On the other hand, for odd values of N, the ground state of 
the spins in the chain is 'nearly' a doublet, and that of the total system must be a doublet. 
Both, the singlet and triplet states of the two impurity spins, can couple to the spins in 
the chain to form a doublet. In fact, for sufficiently large systems in which the gap is much 
smaller than J', the states of the spins in the chain can easily be deformed locally to form 
a singlet with the impurity spins. It seems difficult to obtain reliable limiting values for 
the correlation function from the present sizes of chain when J' is small. 

C. Impurity lattice 

In this subsection, we study systems in which the impurities form a lattice, as shown in 
Fig. 1(c). Each impurity spin is coupled to a spin in the chain with exchange constant J', 
and the impurities do not interact directly with each other. 

We calculate the binding energy per impurity spin which is defined by 

E'^\J\ N) = [Eg{0, N) - Eg{J', N)] /N. (2.3) 

Figure 8 shows the calculated values for E^^^{J' , N). Although the sizes of chain are smaller 
(A?" = 3 ~ 14) than for the single impurity case, the two series for even and odd values of 
N approach each other very closely, and therefore we can extrapolate them accurately to 
AT — > oo. The extrapolated value of E^^^{J') is nearly the same as for A£?i(J'). 

We also calculate the correlation function between one of the impurity spins and the 
nearest spin in the chain. Figure 9 shows the calculated function — < Sq/ ■ Sqc >• For this 
quantity, it is difficult to extrapolate the value to A" — > oo, partly because of small sizes 
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of chain. But even without the extrapolation, we can guess rather accurately the limiting 
value, since the two series for even and odd values of N are already sufficiently close. The 
amplitude of the singlet increases with increasing values of J'. This value is nearly the 
same as for the single-impurity case. When taken together with the result for the binding 
energy, we may state that the local character of each impurity is little influenced by the 
presence of the impurity lattice. 

Inspecting E^^^{J') more closely, we notice that it is slightly larger than AE2^{J')/2, 
as shown in Fig. 5. Forming an impurity lattice leads to an increase in the binding energy. 
Also comparing < So/ ■ Sqc > with that for the single impurity case shown in Fig. 3, we 
notice that the absolute values are slightly larger here. Forming an impurity lattice results 
in an increase in the amplitude of the local singlet. 

Figure 10 shows the correlation function — < So/-Si/ > between a pair of impurity spins 
on nearest-neighbor sites. The dependences on J' and N look complicated, particularly for 
small values of J'. The extrapolated values first increase and then decrease with increasing 
values of J'. For J' > 0.5, the values are very close to those for the two-impurity case 
(Fig. 7), indicating again that forming a lattice affects little the pair correlations between 
impurities. 

In addition to the ground-state properties, we also calculate the energy and the wave 
function of the lowest excited states. The lowest excited states are triplets with momentum 
k = TT. The energy gap A{J',N) is defined by the difference between the energy of the 
lowest excited state and of the ground state. Figure 11 shows the calculated values of 
A{J',N). The gap decreases with increasing N (even and odd), and no lower limit of the 
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gap exists. The extrapolated value A (J') is much smaller than E^^^{J') for small values of 
J' (< 0.5). For J' < 0.3 it is nearly zero. This favors a non-analytic dependence on J' , i.e., 
A(J') oc e-«/-^' with a some numerical constant, as known for Kondo lattices. ■'^■'^ Clearly 
we cannot rule out other dependences on J' or the presence of a critical value of J' below 
which the system is gapless. For larger values of J' (> 0.7), the gap is not much different 
from £'^^(J'). This means that around each impurity a local singlet is well developed, and 
that it is broken up when an excited state is generated. This process requires an energy 
of order S|*(J'). 

Breaking up a local singlet to create excited states becomes apparent in the strong 
coupling limit when J' ^ oo. In that case, each impurity spin forms a perfect singlet with 
the nearest spin in the chain. The low-lying excited states consist of a triplet excitation at 
a particular site with the other sites remaining within the singlet state. They are 3N-fold 
degenerate, having energy J'. The first-order correction to order 1/J' lifts the degeneracy. 
Let \T^{j) > be the state of the triplet to which site j is excited with the magnetic quantum 
number ±1,0). Then the state \ip''ik) > = {^^^'^ J2j^~'''^\T^i3) > is an eigenstate 
with an excitation energy = J' + ^ cos/c.^^ Note that E}~ has a minimum ai k = tt, 
which is consistent with the present numerical results for weak and intermediate couplings. 

Finally, we discuss briefly the changes which occur when the coupling is ferromagnetic 
J' < 0. The ground state is again a singlet. Figures 12 shows the binding energies per 
impurity spin when the impurity number changes. Just like in the case of J' > 0, results 
fall into two groups corresponding to even and odd values of N. If each impurity spin forms 
a perfect triplet with the nearest spin of the chain, the energy gain due to Hfc is || J'|. 
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The extrapolated values are again smaller than ^\J'\, but the reduction ratio is less than 
for J' > 0, as is noticed by comparison with Fig. 8. 

Figure 13 shows the correlation function < Sq/ • Sqc > between one impurity spin and 
the nearest spin in the chain as J' varies. The values first increase rapidly and then move 
gradually toward 0.25 , which is the one of a perfect triplet. By comparison with Fig. 9, the 
local triplet character seems to be established rapidly. Similarly as in Fig. 9, it is difficult 
to extrapolate the value for A'" — > oo. 

The lowest excited states are found to be triplets with momentum k — n. Figure 14 
shows the energy gap A( J', N) between these states. The extrapolated values are very 
small; for —J' < 0.2, the gap is nearly zero. They increase very slowly with increasing 
values of —J'. The gap is much smaller than £'^*(J') given in Fig. 12, indicating that 
the lowest excited states can be formed without destroying the local triplets. In the 
limit J' — > — oo, each impurity spin forms a perfect triplet with the nearest spin in the 
chain. Therefore we can regard each site as being occupied by a spin S = 1 with effective 
exchange coupling constant Jefj = 1/4. For such a system it is known that the lowest 
excited state has momentum k = tt, and a gap which is given according to Haldane by 
A(J') ~ 0.411Jeff = 0.103.^^ Although for J' ~ —1.0 the system is to good approximation 
locally in a triplet state, the gap is found to be still much smaller than the Haldane gap. 

III. CONCLUDING REMARKS 

We have studied systems of magnetic impurities coupled to a spin-| antiferromagnetic 
Heisenberg chain. We have calculated the binding energy per impurity spin and the cor- 
relation function between an impurity spin and the nearest spin in the chain for systems 
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with a single impurity, with two impurities, and with impurities forming a lattice. We have 
shown that the local singlet character increases with increasing antiferromagnetic coupling. 

In the conventional Kondo problem, the binding energy may be written as AEi ~ 
(3 ln2t){pJ')^t+4:te-^/P-^' for weak couplings,^^ where 4t is the band width and p — (27rt) ^ 
is the density of states of the conduction electrons at the Fermi energy.''' The first term is 
the normal part, and the second term is the anomalous part corresponding to the Kondo 
temperature. For t ~ leV and J' ~ 0.2eV, we obtain AEi ~ 0.002eV+4e-^°''eV, i.e., 
the anomalous part is extremely small. On the other hand, in the present Heisenberg 
chain model, the exchange interaction between conduction electrons may be estimated as 
J(= At^ /U) ~ leV when U ~ 4eV and t ~ leV, where U is the Coulomb interaction 
between conduction electrons.^ For J' ~ 0.2eV, we have J' / J = 0.2, and from Fig. 2 
we estimate the binding energy as AEi ~ 0.025eV. This value is at least by one order 
of magnitude larger than that for the conventional Kondo problem. This enhancement 
results from the fact that the disturbance in the chain due to the presence of the impurity 
is much more local for finite U than for t/ = 0. 

We have found that the binding energy per impurity spin as well as various spin cor- 
relation functions are nearly independent of the impurity density. This implies that the 
state which an impurity is forming with its surroundings remains unaffected by the pres- 
ence of other impurities, which indicates its local character. The binding energy is mostly 
determined by this local character. Following the RKKY interaction, we may express the 
indirect interaction between impurities at sites i and j as i^ind — J''^x{h j)^if ' Sj/, where 
xihj) = ~i /o°° < ['^jci'^)i '^ic] > e~^*dt with 5 — > 0+. This interaction was estimated to 
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be small here. Although the interaction energy is small, the spin correlations < Sj/ • Sjf > 
are considerable for small values of J'. They decrease though with increasing values of J'. 
For large J' each impurity spin is forming a singlet with the nearest spin of the chain. 

Clarke, Giamarchi, and Shraiman^^ have recently studied impurity spin models by trans- 
forming them with the help of the bosonization technique to the two-channel Kondo model. 
They have found an anomalous temperature dependence for the susceptibility. Numerical 
studies of such a finite-temperature behavior are left for the future. 
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APPENDIX 

We analyze the case that an impurity is coupled to an anisotropic Heisenberg model by 
using the spin- wave expansion. For an Ising-like anisotropy (a < 1 in Eq.(1.2)), the ground 
state is Neel ordered, and the system is divided into A (up spins) and B (down spins) 
sublattices. In the following, the impurity spin is assumed to be antiferromagnetically 
coupled to a spin in sublattice A, thus pointing downward. We treat the deviation from 
the classical directions by using the Holstein-Primakoff transformation: 

S!, = S-alai, S+ = (S-)^ ^ V2Sa,, (Al) 

S], = -S + b]bj, S+ = {S-)^c^V2Sb], {A2) 
16 



S^f = -S + c^c, S+ = (5o7)t ~ V2Scl {A3) 

Here and bj are boson annihilation operators for spins on sublattices A and B, re- 
spectively, and c is a boson annihilation operator for the impurity spin. Substituting 
Eqs.(Al)-(A3) into Eqs.(1.2) and (1.3), we have 

Ha = -JS'^N + JSJ2 + b]bj + a{a^hj + ajfe])] , (^4) 

<ij> 

Hfc = -J'S^ + J'S(alao + c'c + oqc + ajc'''). (^5) 

Hereafter we measure energies in units of J. 

Since the Hamiltonian is composed only of the terms quadratic in the boson opera- 
tors, we can solve the problem by using a Green's function formalism. ^'^ We consider the 
following functions, 

G{i,i';t) = -i<T{ai{t)al{0))>, D{t) = -i < T{c^ {t)c{0)) >, {A6) 
F{i- t) = -i< T{ai{t)c{0)) >, F{i- t) = -i< T(ct(t)aj(0)) > . (^17) 

Here T is the time-ordering operator. We take He — J'S"^ + J'Sc^c as the unperturbed 
Hamiltonian, and J'S{aQao + qqc + OgC^) ^ ^ perturbation. Then the temporal Fourier 
transform of the unperturbed Green's functions leads to 

k ^ ' 

where 

Efc = 25(l-Q;2cos2A;)V2^ (yllO) 



2Ek J \ 2Ek 
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with — 7r/2 < k < 7r/2. The excitation energy Ej^ has a gap. 

Summing up the diagrams to infinite order, the ones of lowest order shown in Fig. 15, 
we obtain 

G{i, i'; u) = G\i, i'; co) + G\i, 0; j}""^ G\0, i'; a;), {A12) 

1 — ijr [U, U; OJ jTTyOJ ) 

D(.) = + D"M(J'S)^-§^^^^^(.J'S)D»M, (An) 

^' = «°('' "' i-G°(0,'o;.).m '^'"')^°'"'- 

where 

7r(a;) = + {J' Sfo'^iu). {AlQ) 
The energy of the bound state Eb is determined by the relation 

1-G\Q,Q-Eb)ii{Eb)^Q. {An) 

For J' < 2, one of the bound states is positioned in the region —2Sy/l — cP- < Eb < 0, 
and another has in 25' < Eb- Using the Green's functions, we obtain the averages of spin 
operators from the relations, 

< S^f > = -S + ^^lim J D{uj)e-''^*^, {A18) 



27r 



<S^^> = S-i lim [ G(0,0;a;)e-*'^*^, (^19) 
t^o- Jc+ 27r 

<Sl> = S-i f G{i,i; uj)e-''^* ^ . (^20) 

Here C_|_ (C_) indicates that the integration path for du; is taken along the half circle in 
the upper (lower) plane for complex uj. 
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Figure 16 shows the calculated averages of the spin moments for J' = 0.5, S = 1/2. 
The sublattice magnetization < 5'?, > (for i far from the impurity site) decreases with 
increasing values of a, due to the increase of the quantum fluctuation (curve c). The 
spin- wave expansion fails for a close to 1; for a > 0.98, < Sf^. > becomes negative. The 
average of the spin in the chain closest to the impurity site also decreases with increasing 
values of a (curve b). The values are smaller than the sublattice magnetization, since the 
coupling to the impurity spin adds to fluctuations. The absolute value of the average of 
the impurity spin | < Sqj: > \ is reduced from the classical value 0.5, due to transverse spin 
fluctuations caused by the coupling to spins in the chain. They change rather slowly with 
changing a, except when a 1. 

The transverse spin correlation between the impurity spin and the nearest spin in the 
chain are calculated from 

< 5oVo, >= 4 lim / [F(0;a;) + F(0;a;)]e--*^. (^21) 
2 t^o- Jc^ 27r 

Figure 17 shows the values for J' — 0.5, S — 1/2. They are close to the minimum -0.25 
when 1. 
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Figure Captions 

Fig. 1. Sketches of system with (a) a single impurity, (b) two impurities, and (c) impurities 
forming a lattice. Open circles represent spins. Furthermore J and J' represent 
the exchange interactions between spins in the chain and between an impurity spin 
and a spin in the chain, respectively. 

Fig. 2. Binding energy AEi{J';N) for systems with a single impurity. The solid lines 
correspond to even numbers N, from N = A (lowest line) to A?" = 22 (upmost 
line). The dashed lines correspond to odd numbers N, from N = 5 (upmost line) 
to N = 23 (lowest line). The values extrapolated to — > cxd are shown by open 
circles with error bars as indicated. The broken line represents | J', i.e., the energy 
of the perfect singlet. 

Fig. 3. Correlation function — < So/ • Sqc > for systems with a single impurity. Solid 
and dashed lines correspond to even and odd numbers of N, in the same way as 
in Fig. 2. The extrapolated values are shown by open circles with error bars as 
indicated. 

Fig. 4. Binding energy per impurity spin AE2^{J', N)/2 for systems with two impurities 
placed on nearest-neighbor sites. Solid lines correspond to even numbers N, from 
AT = 4 (lowest line) to AT = 22 (upmost line). Dashed lines correspond to odd 
numbers N, from N = 3 (upmost line) to AT = 21 (lowest line). The extrapolated 
values are shown by open circles with error bars. 

Fig. 5. Binding energies per impurity spin, AE2^{J')/2, A£'™(J')/2, and i?^*(J'), in 
comparison with AEi{J'). 
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Fig. 6. Correlation function — < Sq/ • Sqc > between an impurity spin and the spin in 
the chain closest to it, when two impurities are placed on nearest-neighbor sites. 
Solid and dashed lines correspond to even and odd numbers of N, in the same way 
as in Fig. 4. The extrapolated values are shown by open circles with error bars as 
indicated. 

Fig. 7. Correlation function — < So/ • Si/ > for two impurity spins placed on nearest- 
neighbor sites. Solid lines correspond to even numbers N, from N = A (upmost 
line) to A?" = 22 (lowest line). Dashed lines correspond to odd numbers N, from 
N = 3 (lowest line) to A?" = 21 (upmost line). The extrapolated values are shown 
by open circles with error bars as indicated. 

Fig. 8. Binding energy per impurity spin E^^^{J';N) for impurity lattices. Solid lines 
correspond to even numbers N, from N = A (lowest line) to N = lA (upmost line). 
Dashed lines correspond to odd numbers N, from N = 3 (upmost line) to A?" = 13 
(lowest line). The extrapolated values are shown by open circles with error bars 
as indicated. 

Fig. 9. Correlation function — < Sq/ • Sqc > between one of the impurity spins of the 
lattice and the nearest spin in the chain. Solid and dashed lines correspond to 
even and odd numbers of N, in the same way as in Fig. 8. 

Fig. 10. Correlation function — < Sq/ • Si/ > between a pair of impurity spins on nearest- 
neighbor sites, for impurity lattices. Solid lines correspond to even numbers of 
N, from = 4 (upmost line) to A^ = 14 (lowest line). Dashed lines correspond 
to odd numbers N. The extrapolated values are shown by open circles with error 
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bars as indicated. 

Fig. 11. Energy gap A( J', N) for impurity lattices. Solid lines correspond to even numbers 
of N, from = 4 (upmost line) to = 14 (lowest line). Dashed lines correspond 
to odd numbers of N, from N = 3 (upmost line) to A?" = 13 (lowest line). The 
extrapolated values are shown by open circles with error bars as indicated. The 
broken line represents J' — ^, which is the value of A (J') valid for J' » 1. 

Fig. 12. Binding energy per impurity spin E^^^{J';N) with J' < for impurity lattices. 

Solid and dashed lines correspond to even and odd numbers of N, in the same 
way as in Fig. 8. The curve for N=3 is identical to the broken line j\J'\, i.e., the 
energy of the perfect triplet. The extrapolated values are shown by open circles 
with error bars as indicated. 

Fig. 13. Correlation function — < Sq/ • Sqc > between one of the impurity spins of the 
lattice and the nearest spin in the chain, with J' < 0. Solid and dashed lines 
correspond to even and odd numbers of N, in the same way as in Fig. 8. 

Fig. 14. Energy gap A( J', N) with J' < for impurity lattices. Solid and dashed lines 
correspond to even and odd numbers of N, in the same way as in Fig.ll. The 
extrapolated values are shown by open circles with error bars as indicated. The 
broken line represents the Haldane gap for 5' = 1. 

Fig. 15. Diagrams for the Green's functions. Solid lines represent G°(z,z';a;), and broken 
lines represent D^(a;). Crosses represent the interaction J'S at site i = 0. 

Fig. 16. Average of the impurity spin — < Sqj- >, of the spin in the chain closest to it 
< Sq^ >, and of a spin in the chain far away from it < Sf^ >, as a function of a 
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with J' = 0.5. 

Fig. 17. Transverse spin correlation function < Sq^Sq^ > as a function of a with J' = 0.5. 
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